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ABSTRACT 

Context. To investigate and specify the statistical properties of cosmological fields with particular attention to possible non-Gaussian 
features, accurate formulae for the bispectrum and the bispectrum covariance are required. The bispectrum is the lowest-order statistic 
providing an estimate for non-Gaussianities of a distribution, and the bispectrum covariance depicts the errors of the bispectrum 
measurement and their correlation on different scales. Currently, there do exist fitting formulae for the bispectrum and an analytical 
expression for the bispectrum covariance, but the former is not very accurate and the latter contains several intricate terms and only 
one of them can be readily evaluated from the power spectrum of the studied field. Neglecting all higher-order terms results in the 
Gaussian approximation of the bispectrum covariance. 

Aims. We study the range of validity of this Gaussian approximation for two-dimensional non-Gaussian random fields. 
Methods. For this purpose, we simulate Gaussian and non-Gaussian random fields, the latter represented by log-normal fields and 
obtained directly from the former by a simple transformation. From the simulated fields, we calculate the power spectra, the bispectra, 
and the covariance from the sample variance of the bispectra, for different degrees of non-Gaussianity a, which is equivalent to the 
skewness on a given angular scale 9 g . In doing so, we minimize sample variance by selecting the same 'phases' of the Gaussian and 
non-Gaussian fields. 

Results. We find that the Gaussian approximation of the bispectrum covariance provides a good approximation for degrees of non- 
Gaussianity a < 0.5 and a reasonably accurate approximation for a < 1, both on scales > 8# g . For larger values of a, the Gaussian 
approximation clearly breaks down. Using results from cosmic shear simulations, we estimate that the cosmic shear convergence 
fields are described by a < 0.7 at e ~ 4". 

Conclusions. We therefore conclude that the Gaussian approximation for the bispectrum covariance is likely to be applicable in 
ongoing and future cosmic shear studies. 

Key words, cosmology: theory - large-scale structure of Universe - cosmological parameters - methods: numerical - statistical 



1. Introduction 

Huge cosmological surveys are currently ongoing or planned for 
the near future. The large statistical power of these surveys al- 
lows us to study the statistical properties of the cosmological 
fields well beyond the lowest, second-order statistical level. In 
order to employ these higher-order statistics in obtaining infor- 
mation about the underlying cosmology, one needs to be able to 
accurately predict these statistics as a function of the model pa- 
rameters. Moreover, one also needs to predict the expected errors 
on these higher-order measures, i.e., the covariance of higher- 
order statistics. 

Whereas there exist rather accurate fitting formulae for 
the power spectrum of the cosmic d ensity field (e.g., 
iPeacock & Doddsl Il996t ISmifh et all l2003h . this is no longer 
the case for higher-order statistical properties. For example, the 
semi-empiric al fitting formula for the bispectru m of the cosmic 
density field dScoccimarro & Couchma n 2001) is estimated to 
be accurate only at the level of several tens of percent. 

In order to obtain the covariance of the power spectrum and 
the bispectrum, even higher-order statistical properties of the 
field need to be known. The power spectrum covariance depends 
on fourth-order statistical properties, i.e., the trispectrum, which 



is even harder to predict. In principle, one needs to understand 
the sixth-order properties of the field in order to determine the 
covariance of the bispectrum. 

However, depending on the statistical field under considera- 
tion, certain simplifications may apply. In particular, if the den- 
sity field is sufficiently close to a Gaussian field, the higher-order 
terms in the covariance of the power spectrum and the bispec- 
trum may become subdominant, in which case they can be ob- 
tained solely from the power spectrum itself. 

In this paper we study the covariance of the bispectrum for 
two-dimensional random fields. Whereas our prime interest is 
focused on the convergence field of cosmic shear (as a prime 
example of a two-dimensio nal cosmological field) (for a review 
of gravitational lensing, see ISchneiderfl2006l) . our results are not 
restricted to this particular application. 

Third-order cosmic shear statistics is known to pro- 
vide very valuable information, comple mentary to tha t from 
the cosmic shear power spectrum (e.g.. Ijakada & Janl l2004t 
iKilbinger & Schneideijl2005HHuterer et al.ll2006l) . which is con- 
firmed by the first m easurements of third-order cosmic shear 



in previous surveys dBernardeau et aT JgOOl iPen et alj [2003; 



iJarvis et al 1 l2004t ISemboloni et al.l 1201 1). Ongoing and future 
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surveys, such as the Canada-France-Hawaii-Telescope Legacy 
Survey (CFHTLS), the Kilo-Degree Survey (KIDS), the Dark 
Energy Survey (DES) and EUCLID, will measure the third-order 
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cosmic shear signal with very high accuracy. However, estimat- 
ing the statistical power of these surveys, and obtaining realistic 
confidence estimates for model parameters, remains a challenge 
that hinges on our ability to estimate the bispectrum covariance. 

lHul(l2000l) estimated this covariance by analogy to the bis- 
pectrum co variance of the cosmic microwave background. As 
discussed in lJoachimi et al.l (|2009l hereafter JSS), this approach 
has a number of drawbacks that result from the use of spherical 
harmonics in the analysis. Instead, JSS derived a general expres- 
sion for the bispectrum covariance on a finite region in the flat- 
sky limit, which can be separated into four terms. The first term 
contains the contribution of the power spectrum to the bispec- 
trum covariance, and is the only term relevant for Gaussian ran- 
dom fields. The other three terms contain the contributions of the 
bispectrum, the trispectrum and the pentaspectrum. Dropping 
these higher-order terms in the expression is usually termed the 
Gaussian approximation of the bispectrum covariance. If the 
other terms were small for realistic cosmic fields, this would 
tremendously simplify the estimation of confidence regions from 
third-order statistics. At first glance this seems to be invalid as 
the present density field appears to be highly non-Gaussian, es- 
pecially on small scales. Nevertheless, this Gaussian approxima- 
tion has been employed in nearly all studi es of parameter fore- 
casts from third-order cosmic shea r (e.g., iTakada & J ainl 120041 
iHuterer et al.ll2006HShi et al1l2010h . 

In this paper, we investigate the validity of the Gaussian ap- 
proximation, using simulated random fields with different de- 
grees of non-Gaussianity. Specifically, we employ log-normal 
random fields which we characterise by the non- Gaussianity pa- 
ramete r a. As argued in the seminal paper by IColes & Jones! 
(Il99ll) . matter density fields are well described by log-normal 
statistics, maki ng these statistics re levant in a cosmological con- 
text. Indeed, as iKavo et ail ( 1200 ll) have shown, the distribution 
of the density fluctuations in the universe can ra ther well be 
describ ed by a log-normal dist ribution. Moreover, iTaruva et al] 
d2002l) and lHilbertetail (1201 ll) have shown that the same holds 
true for the cosmic shear convergence field, which serves as 
justification for using t hese particular non- Gaussian fields for 
our study. Furthermore, iHilbert et alj d201 ll) employed the log- 
normal ansatz to calculate the covariance of the power spec- 
trum of cosmic shear fields; they showed that this model pro- 
vides fairly accurate results, when compared to ray-tracing sim- 
ulations. We investigate the question to which degree of non- 
Gaussianity a it is justified to use the Gaussian approximation 
for the bispectrum covariance. 

This paper is organised as follows: Sect. 2 reviews the most 
important properties of the analytical formula for the bispectrum 
covariance derived by JSS and the corresponding Gaussian ap- 
proximation, as well as an ansatz to determine the scope of va- 
lidity of the latter applied to non-Gaussian fields. In Sect. 3 the 
primary quantities in the theoretical framework of log-normal 
random fields are introduced, and the method we used to numer- 
ically generate the random fields. The results will be presented, 
analysed and discussed in Sect. 4. Finally, Sect. 5 provides con- 
clusions. 



2. The bispectrum covariance: Analytical formula 
and Gaussian approximation 

In general, the bispectrum of a two-dimensional homogeneous 
and isotropic random field in Fourier space f(f) with zero mean 
is defined by 



and consequently depends on the length t; = \t,\ of the three 
angular frequency vectors l\,t 2 ,t3. It is characterised by the 
fact that the two-dimensional Dirac delta-distribution 6 D (() re- 
stricts the three angular frequency vectors to form a triangle. 
Accordingly, JSS introduced the following bispectrum estima- 
tor for a random field of size Q, 



3 r 

nl 

i=\ Ja r 



^\tx +h + h) , (2) 



so that the bispectrum can be calculated by averaging configura- 
tions over annuli. The mean length of a triangle side, correspond- 
ing to the mean radius of an annulus, is denoted by (, and the 
corresponding area of an annulus is given by Ar(C,) = In^Atj, 
with Aij being the respective bin size. The parameter 



a- 1 [i u t 2 , i 3 ) = - ^2l\t\ + 2~e\e\ + nm 



(3) 



reflects the area of the triangle which is spanned by the three 
wave vectors. JSS showed that the estimator (f2]i is unbiased in 
the limit of small annuli width A^,-. 

JSS derived an analytical expression for the covariance of 
the bispectrum estimator defined in Eq. d2), 

Cav(&tf u h,h), BiUJsJe)) = 7 P(ti)P(h)P(h) 

+ T 3x3 + T 4x2 + T 6 . (4) 

This expression consists of four different terms. The first term is 
proportional to the third power of the power spectrum, with an 
amplitude given by the geometrical factor 



r := 



{2nf 



with the short hand notation 
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and 6xy = 1 if x = y, and zero otherwise. In the case of a 
Gaussian field, this would be the sole contribution to the covari- 
ance; therefore, the approximation of the bispectrum covariance 
by this first term only is called Gaussian approximation in the 
following. The other three terms account for the contribution of 
higher-order spectra to the bispectrum covariance: The second 
term is an integral over the product of bispectra, the third term 
an integral over the product of the trispectrum and the power 
spectrum, and finally the fourth term an integral over the pen- 
taspectrum. 

Due to the complexity of Eq. (0| and the fact that higher- 
order spectra are usually much more difficult to obtain, ei- 
ther analytically or from simulations, usually only the Gaussian 
approximation is used for statistical analyses of cosmological 
fields, despite the fact that these fields are not purely Gaussian 
(see Sect. [1). Thus the question arises to which degree of non- 
Gaussianity of the underlying field it is justified to use the 
Gaussian approximation for computing the bispectrum covari- 
ance. We have addressed the problem by the following approach: 
Applied to a purely Gaussian random field, the Gaussian approx- 
imation should hold per definition, 



(fih) fid) /(/ 3 )> = (2tt) 2 tfg>(/, +e 2 + h) B(i u £ 2 , £3) , (1) Cov G (m , l 2 , ? 3 ), B(U, h, 4)) = 7 ^1)^2)^3) , (7) 
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as will nevertheless be investigated numerically in Sect. 14.3.11 
However, for a non-Gaussian random field, the Gaussian ap- 
proximation might hold to a certain degree of non-Gaussianity a 
(which will be specified in Sect. 13. U as well, 

Cov a (B(£ u { 2 , ? 3 ), B({ 4 , t 5 , 4)) « y Pa{h)P a {h)P a {h) ■ (8) 

In order to assess the range of validity of the Gaussian ap- 
proximation for non-Gaussian random fields, we can compare 
a Gaussian random field with a non-Gaussian one, characterised 
by the degree of non-Gaussianity a, both having the same field 
geometry to assure the same geometrical factor y for both fields. 
The quality of the Gaussian approximation can then be quanti- 
fied, on the one hand, by a parameter g, which is defined by 



Cov a (B(t u € 2 ,li), WuhM) 
Cov G (B(A,4,4),B(A,4,4)) 



JWi)JW2)JW 3 ) 



(9) 



The first term on the right side corresponds to the geometrical 
factor y defined in Eq. ©, the second term to y~ l defined in 
Eq. ©. As long as these to factors give the same result, g{a) « 1 
and the Gaussian approximation provides a reasonable descrip- 
tion of the diagonal part of the bispectrum covariance. On the 
other hand, one needs to check whether the true covariance ma- 
trix is (close to) diagonal, as predicted in the Gaussian approxi- 
mation. 

3. Methods 

3.1. Log-normal random fields 

To estimate the parameter g(a), we consider a particular case of a 
non-Gaussian field, namely a log-normal field. This is motivated 
by its relative simplicity, as well as the fact that the density fluc- 
tuations in cosmological models can reasonably well be approxi- 
mated by a log-normal distribution (see Sect.[T]i. The transforma- 
tion of a two-dimensional homogeneous and isotropic Gaussian 
random field g(x) with zero mean and unit variance to a homo- 
geneous and isotropic two-dimensional log-normal random field 
L(x) is given by 



L(x) = - (e 0<w - c) . 



(10) 



The parameter a controls the deviation of the log-normal field 
from the initial Gaussian field and thus the degree of non- 
Gaussianity of the former. The parameter c and A are chosen as 
such that the log-normal field has zero mean and unit variance, 



dg(x) L(x) p(g(x)) = - (e- /2 - c) = 



=> c - e 



a 2 /2 



(ID 



(L 2 (x))= P 



=> A = 



dg{x) L 2 (x)p(g(x))=A\c 2 -l)=l 
1 



Vc 2 - 1 Ve« 2 - 1 



(12) 



Obviously, c and A depend only on the deviation parameter a so 
that Eq. (TTOb can be reduced in its dependence on a, 
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Fig. 1. Gaussian and log-normal distribution with zero mean 
and unit variance, the latter for different a. The larger a is, the 
stronger the log-normal distribution deviates from the Gaussian. 

The probability density of the Gaussian random field g is a 
Gaussian with zero mean and unit variance, and the one of the 
log-normal random field L can be obtained from the Gaussian 
by the relation p (g) dg = p a (L) dL, yielding 



Pa (L) = 



1 



V2^a(L + A) 



exp 



Ni + i)+ f] 

2a 2 



2 \ 



(14) 



The impact of the parameter a is visualised in Fig.[T] which com- 
pares the Gaussian probability distribution with the log-normal 
probability distribution for different values of a. p a (L) peaks 
at negative L (for a > 0) and shows an extended tail for pos- 
itive L. With increasing a, these features become more pro- 
nounced. Furthermore, the probability distribution vanishes for 
L < L m i„ = -A, as can be easily seen from Eq. (TlOb . 

3.2. Numerical generation of the random fields 

In order to evaluate the parameter g(a), we need to calculate the 
bispectrum covariance of the Gaussian and the log-normal field, 
as well as the power spectra. This has been done by numeri- 
cal simulations: Using repeated realisations of the random fields, 
we can determine their power spectra and bispectra, as well as 
the covariance of the bispectra from the sample variance among 
these different realisations. More specifically, we generated re- 
alisations of a Gaussian random field with a predefined power 
spectrum (see below) which uniquely determines its properties. 
For each realisation of the Gaussian field, we obtained one real- 
isation of a log-normal field applying the transformation ( TTOb to 
every pixel, for each chosen value of a. The Gaussian field power 
spectrum is normalised such that the variance on the pixel scale 
is unity. In this way, the geometries of the Gaussian and log- 
normal fields are the same, and the phases of their realisations 
are the same, which guarantees to minimize sample variance in 
the comparison carried out in Eq. (0. 

A two-dimensional Gaussian random field can be generated 
easiest in Fourier space using the standard procedure, i.e. as- 
suming periodicity and assigning Fourier coefficients with mean 
squared amplitude 



(15) 
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The variance cr^d^l) of the Gaussian random field g(() is 
uniquely determined by the choice of an input power spectrum 
P(\f\) and the choice of a smallest wave number t^, or equiva- 
lent^, the field size SI. It is needed as input for a Gaussian ran- 
dom number generator (see iGalassi et al.l l2009t) . which draws 
separately the real and imaginary parts of ge for each wave 
vector I from a Gaussian distribution with standard deviation 
o"g(KI) k V-P(l^l)- For our purpose, we chose the input power 
spectrum to be of the form 

P{€) oc tit + { y 4 , (16) 

with a normalisation that leads to a unit variance of g, and we put 
t = 1000. The size of the field is fixed to Q = 10° x 10°, leading 
to a smallest wave number of f m [ n ss 36. The field is discretised 
onto N 2 = 1024 2 grid points. 

As the log-normal field is obtained from the corresponding 
Gaussian field by a transformation in real space (see Eq. [Tot , 
the Gaussian field is transformed from Fourier space to real 
space using the Fast F ourier Transform algorithm FFTW (see 
iFrigo & John"sonll2005h . After the log-normal field is calculated 
according to Eq. dTOb . it is then transformed back to Fourier 
space. 

To achieve a good signal-to-noise for the comparison of 
the bispectrum covariance matrix of the Gaussian with those of 
the log-normal random fields, it is crucial to generate a large 
number iV r of realisations of the corresponding field; we chose 
N T = 1000. For each field the power spectrum and the bispec- 
trum are computed for each of these N r realisations and then the 
average over these N r power spectra and bispectra is calculated, 
as well as the variance of the bispectra, which yields an estimate 
of the bispectrum covariance. 

4. Results 

4.1. Power spectrum 

To determine the range of validity of the Gaussian approxima- 
tion according to Eq. (0, the measurement of the power spec- 
trum of the Gaussian and of the log-normal fields is necessary. 
The power spectrum is computed for every realisation of each 



field as the average of the Fourier coefficients \g{\ 2 and \L(\ 2 , re- 
spectively, over the polar angle of t. 

Fig. [2] shows the mean power spectrum obtained by aver- 
aging over the 1000 power spectra for each field in the range 
i ~ 36 - 26000. To determine the bispectrum covariance, the 
power spectrum and bispectrum are measured in bins equidistant 
with At = 1 10 in the range highlighted in Fig.0(^ = 145-2455). 
In the following, we restrict all measurements to the range t 
1 45 - 2455, since for substantially larger I, the computing time 
for the bispectrum becomes quite long. 

The power spectra of the Gaussian and the log-normal ran- 
dom field do not differ noticeably for a = 0.1; in fact, these 
two cannot be distinguished in Fig. [2] However, for a > 0.5, the 
power spectrum of the log-normal field starts to differ from that 
of the Gaussian. Given that the log-normal transformation has 
positive curvature, it increases the amplitude and decreases the 
width of density peaks. This leads to an increase of power on 
small scale, as is clearly seen in Fig. [2] Since the variances of all 
fields are the same on the pixel scale, this increase of power for 
large t must be compensated by a corresponding decrease of the 
power for large scales. 

4.2. Bispectrum 

According to Eq. ([1]), one way to compute the bispectrum of 
a random field is to calculate the 3-point correlator in Fourier 
space and to normalise it afterwards by the corresponding pref- 
actor. We employed such a code, which generates at the begin- 
ning of the computational process a list of all Fourier coefficients 
contained in the different bins defined in Sect. 14. II For every bin, 
the code loops over all combinations of the wave vectors t\ , (2 
and tj, and checks whether the triangle condition is satisfied, i.e., 
^1+^2 + ^3 =0. The wave vectors are further constrained to form 
triangles close to the targeted shape and size, with an allowed 
deviation of At, = 55. 

We first restrict our consideration to equilateral triangles, i.e., 
that all ti are equal. The mean bispectrum of the 1000 realisa- 
tions of the Gaussian and log-normal random fields for different 
values of a is shown in Fig. [3] As expected, the bispectrum of the 
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Fig. 2. The mean power spectrum, averaged over 1000 realisa- 
tions, of the Gaussian random field and the mean log-normal 
power spectrum for different degrees of non-Gaussianity a. The 
error bars are the standard deviation of the 1000 power spectra 
divided by VTOOO. 
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Fig. 3. The mean bispectrum per bin for equilateral triangles, 
averaged over 1000 realisations, of the Gaussian random field 
and several log-normal fields with different values of the non- 
Gaussianity parameter a. The error bars are the standard devia- 
tion of the 1000 power spectra divided by V1000. 
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Fig. 4. The reduced bispectrum B(f)/3P 2 ({) for a Gaussian and 
several log-normal random fields with different a. The error bars 
are calculated by Gaussian error propagation. 
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Fig. 5. The relation between the theoretically calculated geo- 
metrical factor and the numerically computed geometrical fac- 
tor ynum/Ttheo for seven sets of 1000 realisations of a Gaussian 
random field. For clarity the error bars are not plotted. 



Gaussian field is compatible with zero within the statistical ac- 
curacy of our simulations. The size of the error bars increases to- 
wards smaller £ due to the smaller number of modes (or available 
triangles) per (-bin. For the log-normal fields, the bispectrum de- 
viates significantly from zero even for a = 0.1. For small values 
of a the amplitude of the bispectrum increases with a at all val- 
ues of I considered. As seen in Fig. [2] for a < 0.5, the power 
spectra of the log-normal fields are very similar, and therefore 
the increase of the bispectrum with a shows an increasing de- 
gree of non-Gaussianity. For larger values of a, the bispectrum 
no longer behaves monotonically with a at all I, which is related 
to the dependence of the power spectrum on a. To leading order, 
the third-order statistics of a log-normal field is the product of 
an a-dependent factor describing the degree of non-Gaussianity, 
and the square of second-order statistics. Hence, the decrease 
of the bispectrum with a for larger values of a and smaller £ 
reflects the decrease of the power spectrum with a for these £, 
which results from our normalisation condition. For large £, the 
bispectrum monotonically increases with a. 

4.2.1. The reduced bispectrum - Estimating the degree of 
non-Gaussianity 

In order to remove the leading, quadratic dependence of the bis- 
pectrum on the power spectrum, we consider the reduced bispec- 
trum, defined as 

B(( u h, h)[P(£x)P(h) + W<3) + P(h)P((i)Y X , 

which for equilateral triangles reduces to B(£)/3P 2 (£). This re- 
duced bispectrum is shown for different a in Fig.|4] As expected, 
its amplitude grows with increasing a; for a — 1 .5 it is approx- 
imately five times larger than for a = 0.1. Additionally, the re- 
duced bispectrum increases with t, indicating a larger degree of 
non-Gaussianity on smaller scales. 

4.3. Bispectrum covariance 

4.3.1 . Accuracy for the Gaussian case 

Before testing the range of validity of the Gaussian approxima- 
tion for the bispectrum covariance of non-Gaussian fields, it is 
reasonable to first probe the accuracy of the Gaussian approxi- 



mation for a purely Gaussian random field. For this, we calcu- 
late the ratio of the numerically determined bispectrum covari- 
ance and third-order products of the power spectrum, to obtain 
the geometrical factor y num from the simulations, which is then 
compared to the expected value y t heo given by Eq. ©. The value 
of 7num has been determined for seven sets of 1000 realisations 
individually; the ratio y n um/7theo for these seven sets is shown 
in Fig. The ratio fluctuates around one, with higher fluctua- 
tions for smaller wave numbers 6, which is due to the fact that 
the results on large scales are obtained by averaging over fewer 
triangles than on small scales. Only at the smallest value of £ 
considered do we see a significant deviation of this ratio from 
unity. This deviation is not surprising, given the various approx- 
imations that were made for deriving Eq. (0]), the most important 
of which is that the spatial scale is much smaller than the size 
of the field, so that edge effects can be neglected. For £ = 200, 
corresponding to about 1 deg, this assumption is no longer fully 
justified. However, for larger values of £ we conclude that the 
analytical expression for the geometrical factor is accurate. 

4.3.2. Diagonal elements 

The Gaussian approximation for the bispectrum covariance says 
that for equilateral triangles, 

C m = Cov(M,*i,*i), B(i 2 ,£ 2 J 2 j) = yP\li)8- tlh , (17) 

i.e., that the off-diagonal elements of Ce l c 2 are zero, and that the 
diagonal elements follow the Gaussian behaviour. Here, we first 
study the latter, and show the diagonal elements in Fig. [6] For 
a < 0.5, the diagonal elements for the log-normal case show 
no significant difference from the Gaussian case, whereas for 
a > 0.5 the amplitude of the diagonal elements of the log-normal 
covariance matrix decreases in comparison to the Gaussian co- 
variance matrix, reflecting the decrease of the power spectrum 
with a across the I range considered. For large values of a, 
Ca increases again, due to the larger degree of non-Gaussianity, 
where the higher-order terms of Eq. ©, unaccounted for in the 
Gaussian approximation, start dominating the behaviour of the 
covariance matrix. 

To remove the dependence of the covariance from the power 
spectrum, we consider the aberration parameter g(a) defined in 
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Fig. 6. The diagonal elements (l\ = of the Gaussian co- 
variance matrix compared to the diagonal elements of the 
log-normal covariance matrix for different degrees of non- 
Gaussianity a. For clarity, the diagonal elements are multiplied 
by the factor £\, and the error bars are not plotted. 



Eq. {9). This is shown in Fig.|7]for several relatively large values 
of a. For this figure, we study the aberration parameter down to 
roughly eight times the pixel scale. This figure indicates that g is 
of order unity for a < 1.2, but strongly increases for larger values 
of a, indicating a clear breakdown of the Gaussian approxima- 
tion. To study in more detail the transition to this breakdown, we 
show in Fig. [8] the behaviour of g with i for several smaller val- 
ues of a. This figure indicates that for a — 0.5, g ~ 1 within the 
statistical accuracy, and that g as 1 . 1 for a — 1 .0. For a = 1.2, 
the deviations are already several tens of percent, and clearly 
increase towards larger I. Hence we conclude that the diagonal 
elements of the covariance are well described by the Gaussian 
approximation for a < 0.5, on the scales probed, and that the 
deviation from the Gaussian approximation are confined to less 
than 10% for o- < 1. 
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Fig. 7. Aberration parameter g for large degrees of non- 
Gaussianity a. The error bars are computed by Gaussian error 
propagation using ACov = C,y V2/(1000 - 1). 
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Fig. 8. Aberration parameter g for small degrees of non- 
Gaussianity a. 
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ent a. The wave numbers t\ and 62 rise for each panel from the 
left lower corner to the right upper corner. 

4.3.3. Off-diagonal elements 

The validity of the Gaussian approximation furthermore requires 
the vanishing of the off-diagonal elements of the covariance ma- 
trices. We define the correlation coefficient 



(18) 



which by definition is unity for the diagonal elements, and plot 
this in Fig. [9] for different values of a. As expected, the off- 
diagonal elements are compatible with zero for the Gaussian 
case, and are also very small for the log-normal fields with 
a < 1 . For larger values of a, the off-diagonal elements of the 
correlation coefficient become significantly non-zero, clearly in- 
dicating the breakdown of the Gaussian approximation. The am- 
plitude of the off-diagonal elements increases with t\ and 1%, 
i.e., for smaller scales where the non-Gaussianity of the fields is 
stronger. To show this behaviour in more detail, Fig. [TOl shows 
the correlation coefficient r^ 2 as a function of the wave number 
t\ for various values of C%. For a — 1 .0 the off-diagonal elements 
are still almost zero, but increase for a > 1 .0 with large t\ and 

Combining the results from studying the diagonal and non- 
diagonal elements of the bispectrum covariance for equilateral 
triangles, we conclude that the Gaussian approximation pro- 
vides, on the scales considered (down to eight times the pixel 
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scale), an accurate description of the bispectrum covariance for 
degrees of non-Gaussianity corresponding to a < 0.5, and a very 
reasonable approximation even for a < 1. For larger degrees 
of non-Gaussianity, the Gaussian approximation starts to break 
down, and higher-order terms become important for calculating 
the bispectrum covariance. In general, the Gaussian approxima- 
tion works better for large scales than for small scales, as non- 
Gaussianity contributes more on small scales and therefore the 
Gaussian approximation should be worse for this case. 

4.3.4. Degree of non-Gaussianity for cosmic shear fields 

In order to relate our results to future studies of third-order statis- 
tics of cosmic shear fields, we need to obtain an estimate of the 
degree of non-Gaussianity, here parametrised in t erms of a, of 
realist ic sh ear (or convergence) fields. As shown bv lTaruva et al.l 
(120021) and lHilbertet"aTT(l2011h . the cosmic shear convergence 
distribution is rather well described by a log-normal; from their 
results we can estimate characteristic values of a. 

A comparison of th e probability distr ibutions used in 
iTaruva et ail (|2002|) and in lHiibert etafl(l201 ll) with the one used 
for our investigations (see Eq. (TBI ) shows that the probability 
distributions are equivalent, except that we have fixed A to get 
unit variance for the log- normal field. The co nvergence field k of 
ITaruva et all (120021) and iHilbert et al.l (1201 ll) corresponds to the 
log-normal field L used in this paper. Furthermore, the parameter 
kminl or a:o corresponds to the parameter A, which was introduced 
in Sect. 13.11 Finally, their quantity cr is equivalent to the non- 
Gaussianity parameter a. Thus, the degree of non-Gaussianity a 
of the conve rgence field k can di rectl y be estimated from the cr 
computed bv lTaruva etall (120021) and lffilbert et al.l (1201 ll) . 

For a source redshift of z = 1, a smooth ing angle of 8 — 4.0 
and a LCDM-model one can compute from ITaruva et al.l d2002l) 
that the degree of non-Gaussianity of the convergen ce field k cor- 
respon ds to a * 0.437. The result we obtain from lHilbert et al.l 
(120111) is a ss 0.658 for unsmoothed fields. Here the smoothing 
is implicitly determined by the parameters in the cosmological 
simulation, roughly to ~ 3 "5. Evidently, the value of a is not a 
strong function of smoothing scale here. In either case, the non- 
Gaussianity of cosmic shear convergence fields are estimated to 
correspond to a < 1, from which we conclude that the Gaussian 
approximation for the covariance bispectrum is likely to be ap- 
plicable in cosmic shear studies, at least down to angular scales 
of ~ 30". 
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4.4. Non-equilateral triangles 

The results shown so far are based on calculations of the bis- 
pectrum for equilateral triangles, to have just a single parameter. 
However, to get a more complete view on the range of applicabil- 
ity of the Gaussian approximation, the bispectrum must in prin- 
ciple be calculated for all possible triangle shapes. Here, we will 
consider two more shapes of triangles; for the first case, which 
we term 'truncated triangles', we set (\ — lo — t, £3 = 
whereas in the second case, denoted 'sharp triangles', we set 
i\ = €2 = (, €3 = l.5€. For both cases, we calculate the co- 
variance of the bispectrum, as was done in the equilateral case 
before. 

In Fig. QT| we first check the accuracy of the analytical es- 
timate for the Gaussian case, i.e., the ratio y n um/ytheo is shown 
for equilateral triangles (black), truncated triangles (blue) and 
sharp triangles (green). As seen before, the Gaussian approxi- 
mation provides a good description of the bispectrum covariance 
in the case of Gaussian fields. There is a larger scatter around 1 
for truncated triangles, and a smaller scatter for sharp triangles, 
than for equilateral triangles. This is due to the fact that trun- 
cated triangles probe the field properties on larger angular scales 
than equilateral triangles; the opposite is true for sharp triangles. 
Accordingly there are fewer (more) triangles to average over for 
truncated (sharp) triangles than for equilateral ones. 

For a log-normal field with a = 0.5, we plot the aberration 
parameter q in Fig. Q~2]for the three different kinds of triangles. 
As can be seen, there is no significant difference for truncated 
and sharp triangles to the results obtained for equilateral trian- 
gles. This simple example thus indicates that the Gaussian ap- 
proximation for the bispectrum covariance also holds for other 
triangle shapes. We therefore conclude that the Gaussian approx- 
imation can be applied for log-normal fields with a degree of 
non-Gaussianity a < 0.5. In combination with what was said 
above, this implies that for applications to cosmic shear, the 
Gaussian approximation for the bispectrum covariance is suf- 
ficiently accurate for the range of scales studied. 
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Fig. 12. The aberration parameter g(a 
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5. Summary and Conclusion 

In this paper we studied the range of applicability of the 
Gaussian approximation for th e bispectrum covariance, as de- 
rived by iJoachimi et al.l d2009l) . also for non-Gaussian fields. 
We chose log-normal fields as representative for non-Gaussian 
fields, also motivated by the fact that the convergence fields un- 
derlying cosmic shear are well described by a log-normal one. 
Our study is based on numerical simulations of Gaussian and 
log-normal fields with different degrees a of non-Gaussianity. 
Using 1000 realisations of these fields, their bispectrum was cal- 
culated, and the bispectrum covariance was estimated from the 
sample variance. 

In our simulations, a is a non-Gaussianity parameter of the 
transformation dTOb . This formula is applied to the pixel values 
in realisations of Gaussian fields on a grid. For convenience, 
the pixel value variance was chosen to be unity in the Gaussian 
and in the log-normal case. We could as well have chosen the 
pixel variance in the Gaussian and non-Gaussian field to be cr 2 
without changing the essential results of this paper because this 
would have just modified A in Eq. (1 121 . This would have rescaled 
the entire field by a mere constant factor, leaving the skewness 
(Appendix) untouched. 

The following conclusions are drawn for angular scales 
down to approximately eight times the pixel scale. We have 
shown that the bispectrum covariance remains almost diago- 
nal, as long as a < 1. Furthermore, the deviation of the diag- 
onal elements from those predicted by the Gaussian approxi- 
mation are negligible for a < 0.5, and below 10% for a < 1. 
Fr om simulations of c onve rgence fields in weak lensing shown 
in iTaruva et all (120021) and iHilbert et all (1201 ll) . we found that 
realistic values of a are not larger than ~ 0.7 on angular scales of 
about 4"; hence, the Gaussian approximation is expected to yield 
an accurate description of the bispectrum covariance - which is 
fortunate, given the complexity of the full expression of this co- 
variance. 

Whereas the results depend on the assumed statistical prop- 
erties of the non-Gaussian field - here assumed to be log- 
normal - as well as on the shape of the power spectrum, which 
we chose somewhat arbitrarily, we expect that the accuracy of 
the Gaussian approximation for the bispectrum covariance de- 
pends predominantly on the degree of non-Gaussianity (here 



parametrised by a, but we could equivalently use the skewness 
of the field - see Appendix). 
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Appendix A: Skewness of the log-normal field 

Besides the mean and the variance, higher order moments of the 
log-normal random field L can be calculated easily, for example 
the third-order moment 



(L(*) 3 )=A 3 (c 2 -l) 2 ( C 2 + 2), 



(A.l) 



which is of special interest as it is related to the skewness param- 
eter T (3) by: 



-(3) 



(L(*) 3 ) 



(A.2) 



As the skewness depends ultimately only on a, the latter can 
be regarded as to set the skewness of the log-normal field 
L(x), which can also be visualised by the Taylor expansion of 
Eq. (TA~2l i. 



„ 7a 3 29a 5 ^ . 

T (3) = 3a+ + + 0(a 7 ) . 

4 32 V ; 



(A.3) 
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